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Shear viscosity of the gluon plasma in SU(3) YM theory is calculated nonpertur- 
batively, within the stochastic vacuum model. The result for the ratio of the shear 
viscosity to the entropy density, proportional to the squared chromo-magnetic gluon 
condensate and the fifth power of the correlation length of the chromo-magnetic vac- 
uum, falls off with the increase of temperature. At temperatures larger than the 
deconfinement critical temperature by a factor of 2, this fall-off is determined by the 
sixth power of the temperature-dependent strong-coupling constant and yields an 
asymptotic approach to the conjectured lower bound of l/(47r), achievable in = 4 
SYM theory. As a by-product of the calculation, we find a particular form of the two- 
point correlation function of gluonic field strengths, which is the only one consistent 
with the Lorentzian shape of the shear- viscosity spectral function. 



I. INTRODUCTION AND PRELIMINARY ESTIMATES 

The RHIC data on collective expansion dynamics of the quark-gluon plasma can be de- 
scribed by the relativistic hydrodynamics applied to a system with very large initial pressure 
gradients According to these data, particles of different mass are emitted from the quark- 
gluon-plasma fireball with a common fluid velocity, that is a signature of a hydrodynamic 
expansion. Due to a large elliptic flow in noncentral collisions an agreement between 
the experimental data [J and the predictions of relativistic hydrodynamics can only be 
reached when the flow of the QGP-fluid is treated as almost non-viscous. This leads to 
an indication that, in the vicinity of the deconflnement phase transition, the quark-gluon 
plasma produced in the RHIC experiments behaves more like an ideal quantum liquid rather 
than a weakly interacting gas. The mean free path Lmfp of a parton, which traverses such a 
liquid, is much smaller than the inter-particle distance, which is of the order of the inverse 
temperature {3 = 1/T, i.e. [L^^^^/ (3) <ti 1. 
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One can consider for comparison a weakly interacting dilute-gas model of the quark-gluon 
plasma. There ~ (pat)^^ with p and at standing for the particle-number density and 
the Coulomb transport cross section, respectively. Using the standard estimates p ~ 
and at ~ g'^jS^lng'^, where g = g{T) is the perturbative finite-temperature QCD coupling, 
one obtains {L^^^/ P) ~ l/i^g'^lng'^) ^ 1, that strongly contradicts the above-mentioned 
experimental results. One can check Q| that these results could have only been reproduced 
by the dilute-gas model if the perturbative transport cross section, at, were larger by an 
order of magnitude. This inconsistency of the weakly interacting quark-gluon plasma with 
the RHIC data initiated recent calculations of kinetic coefficients in the strongly interacting 
relativistic plasmas. 

Among these coefficients, the one whose values define whether the plasma can be con- 
sidered as weakly or strongly interacting is the shear viscosity rj. It is related to the above 
(-^mfp//?)-ratio via the estimate (rj/s) ~ (Lmfp//?), where s is the entropy density. Accord- 
ing to this relation, the shear- viscosity to the entropy-density ratio, rj/s, becomes smaller 
when the plasma interacts stronger. For instance, for T ~ 200 MeV and the estimated 
typical mean free path Lmfp ~ 0.1 fm, one has (rj/s) ~ 0.1. On the other hand, since the 
mean momentum change Ap of a parton, which propagates through the plasma over the 
distance Lmfp, is of the order of T, the Heisenberg uncertainty principle forbids the ratio 
(Ljafp/P) ~ L^fp ■ (and therefore also rj/s) to be vanishingly small. Up to now, the 
minimal value of l/(47r) ~ 0.08 for the shear- viscosity to the entropy-density ratio has been 
found in A/" = 4 SYM theory j^. It is thus challenging to find other QCD-motivated models 
where this ratio would be that small. Recently, it has been demonstrated that such 
low values of the (r7/s)-ratio can take place even in the perturbative YM plasma, due to 
the bremsstrahlung gg ^ ggg processes. In Ref. j?! it has been argued, though, that the 
perturbatively calculated collisional viscosity is anyhow larger (and therefore subdominant) 
compared to the so-called anomalous viscosity, which is generated by plasma instabilities. 

In this paper, we calculate the (?7/s)-ratio in the gluon plasma of SU(3) YM theory 
nonperturbatively. We obtain the shear viscosity from the Kubo formula, which relates 
the corresponding spectral density p{u) to the Euclidean correlation function of the (1,2 - 
component of the energy-momentum tensor Ti2(x, X4). This method, proposed in Ref. js], 
has been explored in Refs. j^, [l^ with the aim to simulate shear viscosity on the lattice. 
Here we work in the continuum limit and parametrize the Euclidean correlation function of 
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the energy-momentum tensors by means of the stochastic vacuum model [Ufl. This model 
generalizes QCD sum rules by assuming the existence of not only the gluon condensate 
{g'^{F^i,Y) but also of the finite vacuum correlation length This assumption is justified 
by the lattice results on the exponential fall-off at large distances of the two-point correlation 
function of gluonic field strengths HQ, {F^uix)Fl{0)) ~e-^l"l. By virtue of this finding, 
the model manages to quantitatively describe confinement; for instance, the string tension 
reads a oc {g^i.F%Y)- 

Below we will use a finite-temperature generalization of the stochastic vacuum model, 
accessible by implementing the Euclidean periodicity of the X4-coordinate. In the decon- 
finement phase of interest, such a generalization yields for the spatial string tension agiT) 
a formula [l^ pretty similar to its above-quoted vacuum counterpart. This formula reads 
CTsiT) oc /i~^(T) ^5f^(Fj^)^)^, where /i~^(T) is the correlation length of the chromo-magnetic 
vacuum, and ((^^(Fj")^)^ is the chromo-magnetic gluon condensate, which survives the de- 
confinement phase transition. The temperature dependence of the two main ingredients of 
_ju(T) and ((^^(F^^)^)^, can be extracted from the results of the lattice simula- 
16|. 

Since T12 = fi'^Ff^F^^, one a priori expects from the Kubo formula, where the 
(^i2(0)Ti2(x))-correlator enters, that the shear viscosity r] oc (^g'^{F^^'^Yrp- This is a gen- 
eral prediction of the stochastic vacuum model for all the kinetic coefficients, for example 



the mode 
tions 13, 



for the jet quenching parameter [l7|. In fact, according to the Kubo formula, all the kinetic 
coefficients are proportional to the total scattering cross section of the propagating parton, 
which itself is proportional to ((^^(Fj")^)^ in this model H, [l^. Therefore, since the shear 
viscosity has the dimensionality of [mass]^, one can on entirely dimensional grounds expect 
for it the following result: 

At temperatures larger than the temperature of dimensional reduction, T > T^,, jJ^iT) and 
(^g^{F°j)'^y^ are proportional to the corresponding power of the only dimensionful parameter 
present in the YM action at such temperatures, g'^T, i.e. 

/.(T) oc /T, {g\Fi;r)^ oc ig'T)\ 

On the other hand, the entropy density s(T) oc T^, so that 

^ oc /(T) at T> n. 
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In this paper, we quantitatively answer the naturally arising question of whether or not this 
function manages to get below the (47r) "^-threshold at temperatures T < 5Tc, which are 
accessible experimentally and on the lattice. 

The outline of the paper is as follows. In Section II, by assuming an exponential fall-off 
for the two-point correlation function of the energy-momentum tensors (Ti2(0)Ti2(a;)), we 
obtain from the Kubo formula an integral equation for the spectral density p{uj) of rj. In 
Section III, by using for p{uj) a Lorentzian-type ansatz, with the same correlation length 
as that of (Ti2(0)ri2(x)), we explore this equation for the cases of large and small |/c|'s, 
where k is the number of a Matsubara mode. The solution in the large- 1 A; | limit yields the 
range of variation of the numerical parameter a, which enters the initial parametrization 
of (Ti2(0)Ti2(x)). The solution in the small- |/c| limit can only coincide with the large- 1 A; | 
solution for a single value of a from this range. This fixes a completely and makes further 
calculations straightforward. In Section IV, we analytically calculate the shear viscosity rj. 
In Section V, we use this result for rj to numerically find the ratio rj/s. Also in Section V 
we compare the calculated nonperturbative spectral density p{u) with the perturbative one, 
which dominates at large tu's. In Section VI, we discuss the results of the paper, as well as 
possible further developments. In Appendix A, we illustrate the separation of perturbative 
contributions to the Kubo formula from the nonperturbative ones. 



II. KUBO FORMULA FOR THE SPECTRAL DENSITY 



Shear viscosity rj can be defined through the relation 



7] = 71 



dp 



(1) 



1^=0 



where the spectral density p{uj) is a solution to the following integral equation, called Kubo 



formula 
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id 



dojpioj 



cosh \u {xi — I)] 
sinh(cj/3/2) 



rf^x (ri2(0)Ti2(x,X4-/3n)). 



(2) 



Here the sum on the RHS runs over winding modes. We use the Fourier decomposition on 
the LHS of Eq. ^■ 

cosh \oj (x4 — f )] 



sinh(cj/?/2) 



2T.UJ ^ 



fc=— oo 



(3) 
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where Uk = 2TiTk is the fc-th Matsubara frequency. The idea is to have a similar decompo- 
sition also on the RHS of Eq. (I2j). For the implementation of this idea, the following chain 
of equalities is important: 

y ^ ^ ^^^^ = _L y / rfAA"-^e-4i+(-^/^)1+-^-^ = 



k=—oo 



M(3 



20Fr(a) J, 



poo -r<^ 

Jo , 



oo 

e" 

ra=— oo 



Here a > is some parameter, 'T" and "-ft'2-«" stand for the Gamma and the MacDonald 
functions, respectively. We assume that, at T = 0, 

{Tu(0)TM.))=N(a){G'y!^^^l^, (4) 

where {G^) = {g'^^F^yY) and N{a) > is a numerical coefficient, which will be determined. 
Then, in the deconfinement phase (T > Tc) of interest, the above chain of equalities yields 

RHS of Eq. (2) = 2r ■ 7rVr(a)iV(a) (G^)^" J] 2^m2^. ' 

k=-oo + ) 

where {G'^)j, = (^g'^{F°-^'^^^. Note that, for temperatures T > T^,, only the (A; = 0)-term 
in the sum should be considered, since the dimensionally reduced theory is a theory of the 
zeroth Matsubara mode. 

By using also Eq. ([3]), we can now rewrite Eq. ^ in terms of Fourier modes as 

T\/r2a—4: 

"^^^^^J^, = -'^nc)ma) (G% j^^^. (5) 
To solve this equation, we use the parametrization 

p{uj) = ujf{uj), 

where f{uj) is some even function sufficient for the convergence of the integral at large w's. 



Motivated by earlier works [8|, llO|, , we choose it in a Lorentzian-type form 



/H = r- (6) 
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Here C = C(T) > is some function, which will be determined. Apparently, the u- 
integration in Eq. ([5]) converges for any choice of 

a > 0. 

Parametrization of the spectral density through Eq. ([6]) guarantees furthermore that both 



sides of Eq. ([5]) have the same leading large-|fc| behavior. It also implies 



jj that M = M{T) 



is the momentum scale below which perturbation theory breaks down. For this reason, M{T) 
should be of the order of the inverse correlation length of the chromo-magnetic vacuum, /x(T). 
Shear viscosity can finally be obtained by means of Eq. ^ as 

We now solve Eq. ([5]) subsequently in the following cases: |fc| 1 and \k\ ~ 1. 



III. CONTRIBUTIONS TO r/ FROM HIGH AND LOW MATSUBARA MODES 



A. (|A;| > l)-case 



Plugging the Lorentzian-type ansatz ([6]) into Eq. ([5]), we obtain 



LHS of Eq. (5) = C / duj 
Jo 



(o;2 + cu2)(cj2 + M2)'^+l 



u;2" L _L 1 _L a\<^+- 



where x = cu/cj^ and = M/oj^. At \k\ ^ 1, one has 



and the latter integral yields 
LHS of Eq. (5) ^ 



2a 1 c2a 



M < 1, 



4r(a + i) ^^^^^^^^ 



+ 



TT 



(8) 



2 sinfvra' 



We see that, if a < 1, then the leading term of the expansion is oc 



(9) 

i.e. it is k- 



2 sin(7ro) ' 

independent, whereas otherwise the leading term of the expansion becomes /c-dependent. 
For this reason, we restrict ourselves to 



a < 1 
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only. Furthermore, since a^l <^ 1 as well, one can expand in powers of also the RHS of 
Eq. ([5]) to obtain 

RHS of Eq. (5) = 7r22°r(a)Ar(a)^^^^M"[l + O(efc)]- (10) 
Equations ^ and (fTOl) yield the function C: 

C = 7r2"+^r(a)Ar(a) sin(7ra) {G^)^ M^— ^ 
Accordingly, Eq. ((71) yields for the shear viscosity 

= 7r22"+ir(«)Ar(«)sin(7ra)^^^. (11) 
' |fc|»i ^ ^ ^ ^ ^ ^ M5 ^ ^ 

The parametric dependence of this expression on (G^)^. and M is indeed the one following 
from the elementary dimensional analysis made in Introduction. 

B. (|A;| ~ l)-case 

Consider |/c|'s sufficiently small for the inequality 

|cjfc|<M, i.e. 161 >1, (12) 

to hold. Disregarding terms 0{ijj1/M'^) and higher, one has 

C dz r(a) C 



LHS of Eq. (5) 



M^^io (^2 + 1)°+^ 2 r(a + i)M2"' 
where 2; = a;/M, while 

RHS of Eq. (5) ~ T:^2''V{a)N{a)'^^f. 
We obtain from these two equations 

C = 7r3/22-+ir (a + ^) Nia) {G% M^"-^ 
and, according to Eq. (I7|), 

In particular, at T > T*, where only the {k = 0)-mode should be considered, this result 
becomes exact. 
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Figure 1: The ratio 



I |fc|»i 



I |fc|~i 



given by Eq. ([25 



IV. rj FROM THE CORRELATION FUNCTION (Ti2(0)ri2(x)) 



We determine now the parameters N{a) and M, which enter the correlation function ([4]). 
This correlation function reads 

{Tumuix)) = {g'F,^,{0)F,^,{0)FUx)Fl{x)) = 

= {{9%Mm,mU^)Fi{x))) + (^7^F,;(0)F,;(0)> {g'Fl{x)Fl{x)) + 

+ {9'F,^,ml{x)) {g'F^^{^)Fl{x)) + {g'F^,^{Q)Fl{x)) {g'F,^^{0)Fl{x)) , (14) 

where double angular brackets denote a connected (or irreducible) average. We use the 



Gaussian-dominance hypothesis 



which allows one to disregard this connected aver- 



age. For the two-point correlation function of gluonic field strengths we use the standard 
parametrization [nl, 



[g'F;;,{x)Fl{x')) = {6,,6., - 6^ 



d^^Dix 



X 



(15) 



12(iV2 - 1) 

where D{x) is a dimensionless function mediating the confining interaction. In this 
parametrization, we have disregarded a small contribution of non-confining non-perturbative 



interactions 



|26l |. By using Eq. f fTSl) . we obtain for the correlation function f |T4l ): 



{T,2{m2{x)) ^ {g'F^,^{mM {9'Fl{x)Fl{x)) + (/F,;(0)Ft(x)> {g'F^^{Q)Fl{x)) + 
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+ {9'Ft,{^)Fl{x)) {g'F,^,{0)Fl{x)) = ^J^T_^^ D\x). (16) 
The dimensionless function D{x) is usually chosen in the form 

D(a;) = e^'^l^l (17) 

Inserting this expression into the formula for the string tension in the fundamental repre- 
sentation, 

= ^ / d^xD{x), (18) 

one can define the gluon condensate in terms of cTf and the vacuum correlation length /i as 
follows 

(G^> = (19) 

TT 

To obtain for the correlator (Ti2(0)Ti2(a;)) the functional form given by the RHS of 
Eq. ([4]), we modify parametrization (fTTj) to 



D{x)=A{a)^ (2;.|a:|)^-» ' (20) 

where A{a) is a numerical normalization factor. At |x| > /i^^, the new function (l20l l falls off 
with the same exponent as Eq. (fTTll . To obtain the normalization factor A{a), we substitute 
Eq. (l20l) into relation (ITSl) . which holds for any function D{x). Using further expression (fT9l) . 
we obtain 



The correlator (fTBD now reads 



m2(0)T,2(x)) = (G^>^ (22) 



where the function ^(a) is given by Eq. (|2T]) . and we have fixed iVc = 3. Comparing Eq. (l22i) 
with the original definition (111), we conclude that 

N{a) = '^^^ and M = 2/i. 
576 

Equations (fTTl) and (fT3l) yield now contributions to rj from high and low Matsubara modes: 



:2"r(a)sin(7ra) ^ ^ \)'^^ (23) 



|fcl»i 9216 ' ' ' ' /i5(T) 
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Figure 2: Entropy density s{T) in the units of obtained from the lattice values for the pressure 
Plat flfil (courtesy of F. Karsch). 



and 



The ratio of these results, 



A 9216 V 2/ fi^{T) 



(24) 



|fc|>i 



r(«) sin(7ra) 



|fc|~i 



(25) 

v^r(^a + ij 

in the interval < a < 1 of interest is plotted in Fig. [U It equals to unity at a = 1/2, i.e., 
at this value of a, our results for shear viscosity become fc-independent, as they should be. 
This yields the principal analytic result of the present paper: 



7r5/2 [^(1/2) 



(26) 



4608^2 I^HT) 

where ^(1/2) ~ 1.05. Remarkably, ansatz (l6|) at a = 1/2 takes the conventional Lorentzian 
form. In the next Section, we will evaluate the ratio i]/ s numerically. 
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V. NUMERICAL EVALUATION 

Following Ref. {l^, we assume the value Tc = 270 MeV in SU(3) YM theory. We use the 
two-loop running coupling jl^ 

g-\T)=2bo\nj + -^\n(2lnj\, where 5o = h = -' " 



, , , A = 0.104T,, 
3 V 167r2 ' ' 



and Nc = 3 for the case under study. We also assume for /i(T) and for the spatial string 
tension in the fundamental representation, cTf(T), the following parametrizations |l5l. Il7l|: 



/i(T) = /i 



1 at Tc < T < T, 
g^(r)r , T ^ T 



(27) 



cTf(T) = af 



g"(r)T 



92(Td.r.)Td. 



1 at T^<T < T„ 

2 

at T > T*, 



(28) 



131 and af = (0.44 GeV)^. Equation ([T9|), extrapolated to finite 



where /i = 894 MeV 
temperatures, yields for the chromo-magnetic gluon condensate {G^)r 



15 



72 



TX 



af(T)/i^(T) = {G') 



1 at T^<T < T„ 

4 

at T > T*. 



The value of can be estimated from the equation 



0-f(T*) = CTf, 



where crf(T) = \0.566q^(T)T]'^ is the high-temperature parametrization of the fundamental 



spatial string tension 



16l |. Solving this equation numerically, one obtains 



n = 1.28T,. 



The entropy density s = s(T) can be obtained by the formula s = dpia,t/dT, where we 
use for the pressure piat the corresponding lattice values from Ref. [l^- In Fig. [2], we plot 
s(T) in the units of T^. The temperature dependence of the ratio rj/s is determined by the 
function (G^)^ /[/i^(T)s(T)]. One can check numerically that, at T > 2Tc where s/T^ is 
nearly constant, {G% /[fi^{T)s{T)] = 0{g\T)), as was mentioned in Introduction. 

In Fig. m we plot the ratio rj/s, with i] given by Eq. (l26l) . as a function of temperature. 
Also in Fig. m we plot the conjectured lower bound for this ratio, equal to (47r)^^ ~ 0.08, 
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Figure 3: Calculated values of the ratio r//s as a function of temperature. Also shown is the 
conjectured lower bound of (47r)~^ for this quantity, realized in = 4 SYM. 

which is realized in A/" = 4 SYM This bound is indeed not reached by our values, 
although they get very close to it at the highest temperature T = 4.54Tc where the lattice 
data for the pressure (and therefore also for s) are available. 

Furthermore, we compare numerically the obtained nonperturbative spectral density. 



piuj,T) = C{T) 



it) 



where C (T) 



2\2 
T 



uj^ + 4/i2(T) 

with the perturbative one, which at the tree level reads j^, [lo| 

1 



7r\3/2 A\l/2) {G^) 
2/ 576 /l3(r) 



(29) 



207r2 



tanh ^ 



Note that pP^'*(u;,T) ^ only at u; > td{T), where uj{T) ~ 7.5T \m]. For this reason, 
^pert^^ in any case does not affect the calculated rj, which is defined according to Eq. ([T]) 
by the values of the spectral density at u; ^ 0. Figure [4] illustrates the full spectral density 



^fuu = p -|- ppert ^ function of tu/T at T = Tc. At u < ct'(T), p'"^^ is given by the 
obtained result (l29l) . while at u > uj{T) the perturbative part p^^"^ takes it over. Were pP'^''* 
nonvanishing down to u; = 0, it would dominate over p already at u; > 4.85T. That is the 
reason why, hy uo = ijj{T), p^^^^ significantly exceeds p, as one can see from the gap in the 
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values of p at this value of lu. Qualitatively, the same p and the relation between p and 
^pert pgj-gjg^ with the increase of temperature. 

Finally, in Appendix A, we illustrate the correspondence between the splitting of p^^" = 
p + pP^'' and the splitting of {Tu{0)T,2{x)) = {Tu{0)T^2{x)) + {Tu{0)T^2{x)) This 
correspondence allows one to isolate the contribution, which p^'^^*' brings about to the u- 
integral in the full Kubo formula. 



VI. DISCUSSION AND OUTLOOK 

In this paper, we have applied Kubo formula to a nonperturbative calculation of the 
shear viscosity rj in SU(3) YM theory. With the use of the stochastic vacuum model, 
the (Ti2(0)Ti2(x))-correlator entering Kubo formula has been expressed in terms of the 
temperature-dependent chromo-magnetic gluon condensate ((7^(F^")^)^ and the correlation 
length of the chromo-magnetic vacuum p~^(T). As was expected (cf. Introduction), r] turns 
out to be oc /i~^(T) ((^^(Fj")^)^, where the numerical factor is given by Eq. f l26l) . At tem- 
peratures T > 2Tc, the calculated ratio r]/s falls off as 0(^g^{T)), as it should do in the 
dimensionally-reduced theory. Numerically, up to the temperature T = 4.54Tc, where the 
lattice data on bulk thermodynamic quantities are still available, the obtained values of the 
ratio rj/s stay above the conjectured lower bound of (47r)~^, which is reached in A/" = 4 
SYM. 

Formally, our result (l26ll persists even at higher temperatures, being extrapolated to 
which it yields for the (?7/s)-ratio values smaller than (47r)~^. One should, however, re- 
alize that the monotonic fall-off of rj/s with temperature, stemming from the relation 
1] (X p -5(r) {g\F^^^)%, is predefined by our calculational method, which combines Kubo 
formula with the stochastic vacuum model. In fact, all the kinetic coefficients derivable in 
this way should be oc (cf. Introduction). In particular, this is the case for the 

bulk viscosity ( 



which can be obtained from the correlation function 



1 + TTTTT-n TtD^Ix) 



(Here "~" stands for "Gaussian approximation".) On general grounds [23|, one indeed ex- 



pects a monotonic fall-off of the (C/s)-ratio with temperature, as was confirmed by explicit 



calculations 



However, on the same general grounds 



23|, for the {rj/ s)-iatio in question 
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Figure 4: The full spectral density p^'^^^iuj^T) as a function of uj/T at T = T^. At ^ < 7.5, p^"^^ 
is entirely nonperturbative and given by Eq. (f29l) . while at ^ > 7.5 the perturbative part pP^''* is 
dominating. 



one expects the existence of a minimum at temperatures close to [27| and a subsequent 
increase with the further increase of temperature. Indeed, at least at the temperature 
as high as 47ATc, g{T) becomes smaller than unity, and the weakly-interacting dilute-gas 
model of the gluon plasma gradually sets in. As mentioned in Introduction, in the dilute-gas 



model 



25| (rj/s) ~ l/{g'^\ng ^), i.e. this ratio increases with temperature. The stochastic 



vacuum model, on the other hand, being applicable at strong coupling, correctly yields the 
expected fall-off of the (?7/s)-ratio at temperatures < 2Tc, but cannot reproduce its increase 
at much higher temperatures. 

We would also like to emphasize an interesting fact, which has been realized by the end of 
the calculation. We have started with the general a-dependent Lorentzian-type ansatz ([6]) 
for the spectral density p{uj). By using it in the Kubo formula, we have come to the 
conclusion that only for the single value, a = 1/2, this ansate provides the Matsubara-mode 
independence of p{uj). For this value of a, Eq. (I6|) takes the conventional Lorentzian form. 
In this way, also the function D{x) from Eq. ( fTSll is defined unambiguously as 



D{x) =^(1/2) 



K3/2i2fl\x\) 
(2/i|x|)3/2 ' 
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where the value of A{l/2) can be found after Eq. ((26 



Note finally that we have used in our calculation the Gaussian-dominance hypothesis |ll| . 
which allows one to disregard in Eq. (fT4l ) the connected four-point correlation function 
of gluonic field strengths compared to the pairwise products of the two-point correlation 
functions. The same approximation was used in Ref. J2lj| for the calculation of topological 
susceptibility via the four-point correlation function. This approximation can be relaxed in 
the following way. Consider a parametrization for the nonperturbative part of the connected 



22| 



four-point correlation function suggested in Ref. 

{{a'F^u. {x,)f;i^, {x,)f;i^ {x,)f;i^ (x4) > > = 

_ / /^2\'^ j faia2b rasaib f _ \, 

\^ / ] J J V^Mli^lAi3!^3'^/i21'2/i41'4 '^/ill'l/i41'4^M21^2At3i^3/ 



'J J l^Mli^lM4!^4^/i2!^2/i3l'3 ^IJ.lVllJ.2t^2^IJ.3'^'ilJ.il^A) ^ 

, fa2asb ra^a^b / _ \_ 

~rj J lt/ill'l/i2l'2^M3l'3M4i^4 ^Atli^lM3'^3^M2i^2M42^4y' 

g " " 1"/J1M2"1'11'2 "Ml'^2"M2i^l/ l"At3M4"i^3'^4 "/i3!^4 "/i41'3 / ^ 

-f (5 ^ ^(5 ^ ^ [5 ^^pi_^5y-^^yj^ '^/^lZ/4'^/i4Vl) ('^/i2/i3'^i^21'3 ^ fi2l^3^ fiSI^2) J'-^l^l; ■ ■ ■ ; ^g) 1 (30) 

where Zi = Xi — X2, Z2 = Xi — X3,..., Zq = X3 — X4 are relative coordinates. For the connected 
correlation function entering Eq. (fT4ll this parametrization yields 

9%^,{0)F,-^{0)Fl{x)Fl{x))) =2{N'^ -1) (n^-I) {G')' D{0,x,x,x,x,0). 



Similarly to Eq. (I20ll . for the function D too one can have a parametrization compatible 
with Eq. (11): 

D{0, X, X, X, X, 0) = B{a) 



(4^|x|)2-° 

The normalization factor B{a) should now be determined simultaneously with the normal- 
ization factor A{a) from a system of equations for two observables, which both depend on 
the functions D and D. Natural observables of this kind are the string tension and the 
topological susceptibility. The contribution of the function D to the string tension has al- 



ready been evaluated in Ref. 



22| . Further analysis of the outlined extension of the Gaussian 



approximation is, however, not the purpose of the present paper. 
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Appendix A. Matching perturbative contributions in the Kubo formula. 

For this Appendix, we promote Eq. ^ to the full Kubo formula, i.e. replace p by p^"^^ = 
p + pP-* and (Ti2(0)Ti2(x)) by (Ti2(0)Ti2(x))f^ii = {Tu{0)T,2{x)) + {T.^muix)) On 
the LHS of such a full Kubo formula, consider the integral containing Ppert- To facilitate the 
^-integration, we approximate the perturbative part of the spectral density by the function 
ppert _ jY^4 (Jq^j^ to u = and determine the coefficient A^. The integral emerging on the 
LHS of the Kubo formula then reads 



poo ^ 

/ dcuu;^ [cosh(u;a;4) coth(a;/3/2) — sinh(u;a;4)] = 
Jo 



1 1 

+ 



{(3n + Xif [(3n+{(3-Xi)f 
1 



(Al) 



n=l 



{(3n + ^4)^ {[3n — x^Y 
To obtain the last expression, we have extracted the {n = 0)-term from the sum ^ (i3n+x4)^ 

n=0 

00 

and shifted n by 1 in the sum Yl [^n+{^-x4)]^ ■ Note that X4 G [0,/5], and the obtained 

n=0 

expression has singularities at X4 = and X4 = P [in the {n = l)-term]. These singularities 
are identical, since cosh [u (x4 — |)] = cosh(c<j/9/2) both at X4 = and X4 = f3. 

We will demonstrate now that this expression corresponds to the contribution, which 
(^i2(0)Ti2(x))pgj.^^ brings about to the RHS of the full Kubo formula. The UV-finite part of 
this perturbative correlation function can be written as 

AiN^ - ly 



(Ti2(0)Ti2(x)) 



pert ^8 



where the value of the numerical constant A depends on the regularization scheme applied. 
Thus, the sum emerging on the RHS of the Kubo formula reads 



E 



[x^ + {(3n + X4 
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Doing the integration over (fx first, we have 

„2 



dx-. 



X 



TT 



[x^ + {pn + x^y]^ 32|/3n + X4|5' 
The part of the sum with negative winding modes reads 

-1 -, -1 



1 



y L 



1 



\l3n + xA^ ^ (—dn — Xi)^ ^(Bn — XiY 

^ n=— CO ^ ' n=l ^ ' 

Therefore, the integration over c?|x| and the summation over winding modes yield 

1 T^{^ ^ r 1 1 



dxx"^ 



[x2 + (/5n + a;4)2]^ 



(/3n + x^)^ {l3n — X4) 



{A.2) 



^4 n=l 

The coincidence of curly brackets in Eqs. (A.l) and (A. 2) proves that the ansatz p^^^^ = Nu^ 
captures the contribution of (Ti2(0)Ti2(x))p^j.^ correctly, with the corresponding normaliza- 
tion factor being N{T) = 7r^A{T){N^ - 1)^^(T)/192. 

In this way, one isolates in the full Kubo formula simultaneously the perturbative cj^-part 
of p^"" and the perturbative contribution to (Ti2(0)Ti2(a;))£^[[. The remaining nonperturba- 
tive parts of p^"" and (T'i2(0)Ti2(x))£^[[ are related to each other by means of Eq. (I2|). 
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